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1. Introduction 

The cost of the computation of the equation of state (EoS) increases very quickly with in- 
creasing temporal extent A^, and aspect ratio Ns/Nt fixed. Since the lattice spacing for small Nt 
- a = \/{TNt) - can be quite large, especially at low temperatures, it is very important to use 
improved actions with small discretization errors. The MILC collaboration has a longstanding 
program of full 2+1 flavor simulations using the asqtad quark action [|l|] combined with a one- 
loop Symanzik improved gauge action [Q]. This includes a computation of the EoS at zero baryon 
chemical potential and hence vanishing baryon density [Q]. Here we will present the results of an 
extension of this computation to better approximate conditions in heavy ion colUsion experiments, 
namely the inclusion of a (small) chemical potential. 

To avoid the notorious sign problem - the fermion determinant becomes complex with a 
nonzero chemical potential making straight forward Monte Carlo simulations impossible - we use 
the Taylor expansion method iQ] which requires only simulations at zero chemical potential and 
only on the finite temperature ensembles. 



2. Technicalities 



In the Taylor expansion method |Q] one expands, for example, the pressure as 
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Here Z is the partition function, and /i/ /, are the chemical potentials for the Ught and heavy quarks, 
respectively. The expansion coefficients are evaluated at zero chemical potential /X/ /, = 0. Due to 
CP symmetry the terms in the series with n + m odd vanish. The nonzero coefficients are 
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where now the /X; /, are the chemical potentials in lattice units. Knowledge of the Cnm{T) also allows 



for the computation of quark number densities and susceptibilities from eq. ([2. l]), since 
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Similarly, the interaction measure is expanded as 
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where again only terms with n + m even are nonzero and 
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For explicit expressions of the Taylor expansion coefficients Cnm{T) and bnm{T) we refer the reader 
to il. 

To determine the Taylor expansion coefficients Cnm{T) and bnm{T) in numerical simulations, 
we need to calculate traces of derivatives of the asqtad fermion matrix such as 
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and products of such traces. The number of such terms increases fast with increasing order n + m. 
Up to sixth order 40 such terms need to be computed. These traces are estimated on the ensembles 
of lattices along a trajectory of constant physics using 200 random sources in the region of the phase 
transition/crossover and 100 sources outside that region. With these numbers of random sources 
the noise in the Taylor expansion coefficients is dominated by configuration-to-configuration fluc- 
tuations. Increasing the number of random sources would thus not decrease our statistical errors 
substantially. 

3. First numerical results 
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Figure 1: Some nonvanishing coefficients in the Taylor expansion of the pressure as function of the temper- 
ature along the trajectory of constant physics. 



We have computed the Taylor expansion coefficients up to sixth order on the lattice ensembles 
along one trajectory of constant physics where we had computed the EoS at vanishing chemical 
potential before ||3|]. The trajectory is given by the choices Nt = 4 and niud ~ O.lm^, where is 
tuned to the physical strange quark mass to within about 20%. 
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Figure |l] shows our result for some of the coefficients for the contribution to the pressure 
difference A/j = 7^ 0) — p{pii.h = 0), see eq. (2J_). Note that the coefficients quickly reach 

the continuum Stefan-Boltzmann limit above Tc. Also note that the mixed coefficients with both 
n,m^O are quite small. Some coefficients contributing to the difference for the interaction measure 
AI = 7^ 0) —/(/X/./, = 0) are show in figure]^ Since the pressure can be obtained by integration 
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Figure 2: Some nonvanishing Taylor coefficients in the expansion of the interaction measure. 

of the interaction measure along the trajectory of constant physics, a^p = Jj.f^pa"^I{a')dlna' , the 
coefficients c„,„{T) can be obtained by integrating bnm{T), giving a consistency check. This is 
illustrated in figure ^. 
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Figure 3: Comparing the two methods to compute 
C2o(r), directly or by integrating b2o{T). 
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Figure 4: The induced strange quark number den- 
sity ns/T^ when turning on only p-i/T and keeping 
fih/T = 0. 



4 



The equation of state with nonzero chemical potential for 2+1 flavors 



U. M. Heller 



Turning on /i//r induces a small negative strange quark number density ris/T^ because some 
of the c„i(r) are nonvanishing, as shown in figure To keep ris = 0, as is the case in heavy ion 
experiments, we need to compensate by a small tuned p-h/T. We show the pressure and energy 
density contribution due to the nonzero chemical potential, with tuned = 0, in Figure ||. 
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Figure 5: Pressure, Ap, and energy density, Ae, contribution for the m„^ « QAnis, Nt ~ A trajectory to 



4. The isentropic equation of state 

Heavy ion collision experiments produce matter that, after thermalization, is expected to ex- 
pand, at fixed baryon number, without further entropy generation, i.e. isentropically - with constant 
s/ub. At the AGS, SPS and RHIC s/ub is approximately 30, 45 and 300 [^, respectively. To ac- 
count for this situation, we numerically determine the jXi and /i/,, as function of T, by solving 



riB 



(4.1) 



for C = 30, 40, and 300, within our statistical errors. With the determined jU/ and jXh we can then 
compute the isentropic equation of state, shown in figures ^ and ^ (left). For comparison, we also 
include the results with jLt/ = /X/, = 0, i.e. s/hb = 
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Figure 6: Isentropic versions of the interaction measure (left) and pressure (right). 



In figure ^ (right) we show, as a further example, the isentropic light-light quark number 
susceptibility, Xuu- We note that Xuu does not develop a peak structure on any of the isentropic 
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Figure 7: The isentropic energy density (left) and light-light quark number susceptibiUty (right). 



trajectories, as would be expected near a phase transition point. Tiierefore all the mentioned heavy 
ion experiments take place far away from a possible critical (end) point in the jJ. - T plane of the 
phase diagram. 



5. Conclusions 

We have extended the computation of the QCD equation of state for 2 + 1 flavors along a 
trajectory of constant physics with m„d /m^ f» 0. 1 on lattice ensembles with Nt =4 to small nonzero 
chemical potential with the Taylor expansion method up to sixth order in the chemical potential. We 
tuned the strange quark chemical potential to keep the strange quark density vanishing at different 
values of the light quark chemical potential. 

We have also determined the isentropic EoS and quark number susceptibilities for values of the 
ratio s/nB relevant for heavy ion collision experiments. We found no signs of a possible phase tran- 
sition along any of the considered isentropic trajectories. Qualitatively our results are in agreement 
with previous two-flavor studies. 
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